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Abstract: We present analytic expressions of all integrals required to complete the explicit eval- 
uation of the real- virtual integrated counterterms needed to define a recently proposed subtraction 
scheme for jet cross sections at next-to-next-to-leading order in QCD. We use the Mellin-Barnes 
representation of these integrals in 4 — 2e dimensions to obtain the coefficients of their Laurent 
expansions around e = 0. These coefficients are given by linear combinations of multidimensional 
Mellin-Barnes integrals. We compute the coefficients of such expansions in e both numerically 
and analytically by complex integration over the Mellin-Barnes contours. 
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1. Introduction 



Precision predictions in perturbative Quantum Chromodynamics (QCD) at colliders demand cal- 
culating physical observables beyond leading order (LO) accuracy and, in the traditional approach 
to higher order predictions with fully differential kinematics, real and virtual corrections are sep- 
arately evaluated. Integration over the phase space then requires a consistent treatment of the 
infrared singularities before any numerical computation may be performed. At next-to-leading 
order (NLO), infrared divergences can be handled using a subtraction scheme, which exploits 
the universal structure of the kinematical singularities of QCD matrix elements. The necessary 
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(process-independent) counterterms regularize the virtual corrections at one loop and the real 
emission phase space integrals simultaneously [1] . 

At next-to-next-to-leading order (NNLO), the calculation of the radiative corrections to fully 
differential cross sections is a challenging problem and various extensions of the subtraction 
method at NNLO have been proposed, see e.g. Refs. [2-5]. Currently, the available results 
for electron-positron annihilation at NNLO include total rates [6-8] and event shapes [9, 10] for 
the process e + e~ — > 3 jets and are all based on the antenna subtraction method [11-13]. On the 
other hand for colorless final states, such as vector boson or Higgs boson production at hadron 
colliders dedicated subtraction schemes at NNLO [14, 15] have been applied. The infrared struc- 
ture of scattering processes with three or more colored partons is involved if calculated at NNLO 
with the antenna subtraction method [16] - a fact which has motivated the formulation of al- 
ternative subtraction schemes. In particular, Refs. [17-19] introduce a scheme for computing 
NNLO corrections to QCD jet cross sections for processes without colored partons in the initial 
state and an arbitrary number of massless particles (colored or colorless) in the final state. Very 
recently, following the steps of Ref. [17], this subtraction scheme has been extended to cross sec- 
tions for hadron-initiated processes [20], although yet to NLO accuracy only, but in a way which 
is NNLO-compatible. 

Any subtraction scheme is of practical utility only after the counterterms for the regularization 
of the real emissions are integrated over the phase space of the unresolved partons. In the scheme 
of Refs. [17-19] these counterterms are universal and, therefore can be computed once and for 
all. Their knowledge is necessary to regularize the infrared divergences appearing in the virtual 
corrections. Some of the integrals needed explicitly in the so-called real-virtual counterterms of 
this scheme have been calculated in Refs. [21,22]. In the present paper we complete this task by 
computing all integrals needed for the the real-virtual counterterms in the subtraction scheme of 
Refs. [17-19] by means of Mellin-Barnes (MB) representations. The use of MB integrals when 
dealing with Feynman integral calculus has proved powerful in the last years. MB integrals were 
first applied to Feynman integrals in Refs. [23, 24] and pioneering work has been performed since 
then in Refs. [25-27] (see also Ref. [28] and references therein for many other examples). For a 
given integral the MB representation replaces the power of a sum in the integrand by a product 
of the individual terms of the sum raised to some other powers. This leads then to integrals over 
certain complex contours of T-functions. As a crucial point it is then very convenient with this 
MB representation to resolve all singularities in the limit e = within dimensional regularization, 
d = 4 — 2e. In this paper, we adapt the MB method to derive analytic expressions for all integrals 
appearing in the real- virtual counterterms of Refs. [17-19]. 

Let us briefly discuss the merits of the analytic approach for the computation of the inte- 
grated subtraction terms. First of all, in a higher-order computation, the e poles of the integrated 
subtraction terms need to cancel the corresponding e poles coming from the loop matrix elements 
in the virtual corrections. The cancellation of these poles can be demonstrated most convincingly 
once the pole structure of the integrated subtraction terms is exhibited analytically. Second, in 
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terms of speed and precision of the evaluation, analytic results are very fast and very accurate 
compared to numerical ones. Moreover, they demonstrate that the final result consists of smooth 
functions only. Nevertheless also the numerical evaluation of the integrated counterterms has its 
utility, because it serves as an independent check. Then, there are indeed some cases, where it 
is very difficult to find the analytic computation of the multi-dimensional MB integral and only 
the complex numerical integration can be carried out. In these cases, however, the method of 
MB integrals provides a fast and reliable way to obtain the final results with small numerical 
uncertainties. From a practical point of view, the combination of both, analytic and numerical 
evaluations of all MB integrals implies that the final results for the integrated real- virtual coun- 
terterms can be conveniently given e.g. in the form of interpolating tables which can be computed 
once and for all. This suffices for any practical application, because in an actual computation 
the relative uncertainty associated with the numerical phase space integrations is generally much 
greater than that of the integrated subtraction terms. 

The outline of the paper is the following. In Sect. |2] we briefly review the phase space integrals 
of the real-virtual corrections at NNLO and we define the integrals of the subtraction terms that 
we will consider in this paper. In Sect.^ we present a brief explanation of the method of MB 
representations. We outline the steps of our calculation and we also discuss explicitly an example 
to display the typical structure of the integrals we are interested in. In Sect.|4] we complete the 
analytic evaluation of all integrals needed for integrated collinear counterterms. Next, in Sects. |5|- 
we compute also all different types of the nested integrals. Finally in Sect. |8| we present the 
conclusions of this work. 



2. Integrals needed for the integrated subtraction terms 

The subtraction method developed in Refs. [18, 19] relies on the universal soft and collinear 
factorization properties of QCD squared matrix elements. Once the subtraction scheme is defined, 
one has to integrate the subtraction terms over the factorized phase space of the unresolved 
parton(s). This is the content of the present work (see also Ref. [21]). 

There are two crucial elements in the formulation of a subtraction scheme beyond NLO. 
Firstly, the factorization formulae should disentangle the overlaps in soft-singular factors and 
collinear singularities in order to avoid multiple subtractions and a simple solution to this problem 
has been given in Ref. [29]. Secondly, because the factorization formulae are valid only in the 
strict soft and collinear limits, they have to be extended to the whole phase space. Typically, 
this requires a mapping of the original n momenta {p} n = {pi, . . ■ ,p n } in an ra-parton matrix 
element at any order in perturbation theory to m momenta {p} m = {pi, ■ . . ,p m } in such a way, 
that momentum conservation is preserved. Here m denotes the number of hard partons and n — m 
is the number of unresolved ones. 
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Figure 1: Graphical representations of the momentum mappings and the implied phase space factorization: 
collinear (left) and soft momentum mapping (right). 

The original n-particle phase space of total momentum Q reads 

" id ( n \ 

dMPi,---,Pn-,Q) = Uj^p [ 5 + (pt)(2n) d sW \Q-J2Pij , (2-1) 
and, for a given mapping, one obtains the phase-space factorization as 

#»({?}„; Q) = d<M{p} m ; Q) [dp n - m; m(Mn- m ; Q)] , (2.2) 

which was first introduced in Ref. [1] in the context of computing QCD corrections at NLO. In 
this paper we are concerned with the integrals of the singly-unresolved counterterms (i.e. the case 
m = 1), which imply two types of mappings: 

Mn — ^ {P>»-1 = {Pl> • • • >Pir, • • • ,Pn} , (2-3) 
Wn^Mil = {Pl,---,Pn}- (2.4) 



In the collinear momentum mapping — ^> in Eq. (2.3) the momenta p^ and pr are replaced by a 



single momentum p^, and all other momenta are rescaled, while for soft-type subtractions, — ^ in 
Eq. ( |2,4| ) the momentum , that may become soft, is missing from the set, and all other momenta 
are rescaled and transformed by a proper Lorentz transformation. Both momentum mappings and 
the corresponding factorization of the phase-space measure are represented graphically in Fig. |l|, 
where the symbol ig> stands for the convolution as implied by Eq. (|2.2| ). The integration of the 
singly-unresolved subtraction terms requires three basic types of integrals over the corresponding 
factorized phase space, as well as iterations of these (nested integrals are denoted by a *). All 
necessary integrals were derived in Refs. [21,22]. 

2.1 Basic integrals 

The three basic integrals are those used in the collinear, soft and soft-collinear subtraction coun- 
terterms. The collinear integrals have the general form 

l(x; e, a , <k; «, k, 5, g { j>) =x da a -l-( 1 +*)« (1 - a) 2 *" 1 [a + (1 - aja?] -1- ^ 6 

Jo 
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If I x^J 





9A 


1 


Tl 




(1 - 







(l-z) ± ' 2 F 1 (±e,±e,l±e,z) 


±1 




2 F 1 (±e,±e,l±e,l-z) 



Table 1: The values of S and gj{z r ) for which Eq. ( |2.5| ) needs to be evaluated. 



These integrals need to be known as a function of x £ [0, 1] (see Refs. [17-19] for its kinematic 
definition) in a Laurent-expansion in e for k = —1,0,1,2. The necessary values of 5 and the 
expressions for the functions are given in Tab. |l[ Here k = 0, 1 for the first row and k = 1 
for all other cases. Analytic expressions for all cases corresponding to the first two rows of Tab. |l] 
were derived in Ref. [21] and contain the first five terms in the e-expansion. We compute all 
cases anew and present our results explicitly in Sect.|||. The other parameters ao £ (0, 1] and 



do in Eq. fl2.5|) will be specified in Sect.y. Our analytic results for these integrals include all the 
coefficients of the poles in e and the first three terms in the e-expansion. 



Next, the soft subtractions require the integral 



r 2 (i 



ik.C 



2vrr(l - 2e) 



^ 1+K£ ' 1+K£ )(cosx) 



(2.6) 



Jo 



as a function of Y.^ g 6 [0, 1] (see Refs. [17-19] for its kinematic definition) in a Laurent expansion 
around e = 0, where Q( l ' k \cosx) denotes the angular integral in d-dimensions 



n {i ' k \cosx) = d(costf) (sintf)~ 2e / d(cos tp) {simp) 



-l-2e 



x (1 — cos i?) l (l — cos x cost? — sinxsin#cos<£>)~ 



with 



cosx 



1 - 2 Kr 
IK.t 



Finally, the soft-collinear subtractions lead to the integral 

ryo rl 
JC(e,y ,d' -K) = 2 J dy y^ 2+ ^(l - y)^ 1 j d(cos0) (sintf)" 



2c 



1 



2(1 -y) 

y(l — cos??) 



r 2 (i 



2vrr(l - 2e) 



d(cos <£>) (sin 99) 



-l-2e 



1 



(2.7) 



(2i 



(2.9) 



which does not depend on kinematical variables. The integrals J and /C in Eqs. (|2.6|) and ( |2.9|) have 
been computed in Ref. [21] for all relevant parameters in yo, <i and k. We have evaluated these 
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soft and soft-collinear integrals, too, and we have checked that the two results agree numerically. 
We do not deal with the cases J and fC in this paper. 



2.2 Nested integrals 

In a NNLO computation, also iterations of the above integrals appear. In this paper we complete 
the list of nested integrals necessary for the integrated real-virtual counterterms, in particular we 
cover all cases which have not been addressed in Ref. [21]. 

Of the nested integrals, which we generally denote by a star *, three are collinear integrals 
with one of the basic types in its argument, 



rao rx 

T*Ti(x] e, a , d ; k, I) = x / da dv oT 1 ~ t (1 - a) 2d °~ l [a + (1 - a)x] 

Jo Jo 

a + (1 — a)xv 



-l-e 



x [v(l — v)]~ 



2a + (1 — a)x 



(a + (1 — a)x(l — v) 
x — o — m \ ! e ' a o, d ; 0, 1, 0, 1 
2a + (1 — a)x 



(2.10) 



I*I r (x;e,ao,do;k,l)=x / da dva" 1 " 6 (1 — a) 0-1 [a + (1 — a)x\ 

Jo Jo 



l-€ 



x [v(l — v)]~ 



a + (1 — a)xv 
2a + (1 — a)x 



'( a + (1 — a)xv \ , . 

1 \ x ^T^n V; e > a °,do;0,Z,0,l , 2.11 

\ 2a + (1 — a)x J 



which we need for k, I = —1, 0, 1, 2, and 



2*J(x;e,ao,do,yo,d ;k) =x da dv a" 1 " 6 (1 — a) 2d °~ 1 [a + (1 — a)x] 

Jo Jo 



-l-e 



x[v(l — v)] 



a + (1 — a) xv 
2a + (1 — a)x 

a(a + (1 — a)x)(2a + (1 — a)x) 2 



(2.12) 



(a + (1 — a)xv){a + (1 — a)x(l — v)) x- 



yo,d' ,o 



0' U I ! 



for A; = —1, 0, 1, 2. Both, X*X and X*^ are needed as a function of x G [0, 1] in an e-expansion with 
X and J given in Eqs. Q2.5| ) and (|2.6|), respectively. A discussion about the choice of the relevant 
parameters ao, do, yo and d' is given at the end of Sect. |3| and details of the computation are also 
given in Sect. B 



Three other iterated integrals are defined as soft integrals with other soft integrals appearing 
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in the argument, 



r 2 (l- e ) f 1 
J*Jik(Y- kQ ;e,y ,d' ) = -8Y rk Q _ — J ^ 



d(cost?) (sin??) 



-2c 



J d(cosip){smip)- 1 ~ 2e (l-cos'&)~ 1 

\ dy g/~ 1-2e (l — y) °[2 — (1 + cos x) cost? — sin x sin t? costp] ~ x 
Jo 



and 



2h ^ /-l 

d(cost?) (sint?)~ 2e 



i 



x / d(cos ip) (sin tp) 1 2e (l-cost?) 1 



l 



(2.13) 



x ^/ 4 ( 1 -^)%,Q . £ yo (/ , (J 

\ [2 — y(l + cos 1 !?)] [2 — y(l + cos x cost? + sin x sin t? cos cp)] ' ' ' °' 

M*rV*& ^ »>> d 'o) = ~ 8 %Q 2^(1 te) j_ x d(c ° S (sin??r2£ 

x f d(cos (^(sincp) -1-26 (1 - cost?) -1 
J-i (2.14) 

x [2 — (1 + cosx) cost? — sinxsin^cos^] -1 



(2.15) 



rvo 

x / dyy _1 ~ 2e (l — y) [2 — (1 + cosx) cost? — sinxsint?cos( / 9] -1 
Jo 

( (1 — cosx cost? — sin x sin i? cos cp) ,/ 

x J\ 7, T\ | — q -r;e,yo,a ,0 

\ 2 — y(l + cos x cos if + sin x sm v cos 0) 

with J given in Eq. ( |2.6[ ). The three integrals in Eqs. fl2.13| )-( p?.15| ) need to be calculated for 
Y-q. q £ [0,1] as expansion in e. Explicit results and details of the computation (and values for 
the parameters yo and d' ) for these integrals are presented in Sects. |^ and ^. 

The final case is when the soft integral appears in the argument of a soft-collinear one, 

/C^(6,y ,rfo) = 2 2 ^ ( ( 1 l ~ £ 2 ) e) J d(cc*0) (sintf)- 2e 

/i rvo , 

d(cosvj)(sin^)- 1 - 2e J dyy- l - 2e (l-y) d °- 1 (2.16) 

2-t/(l + cost?) / 1-cost? \ 
1 — cost? \2 — y(l + cost?) / 
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which is again independent of the kinematics, i.e. the coefficients of the expansion in e are pure 
numbers. Details of the computation and the parameters yo and d' are given in Sects. |3| and [?]. 



3. The method of Mellin-Barnes representations 



In this Section we briefly review the essential steps in the derivation of MB representations for 
the integrals of Sects. EO and 2.2. The starting point is the well known basic formula, 



1 



1 



(a + bf I» 



q+ioa j 

— a- u - z b z T(u + z)T(-z) 

q—ioa 



(3.1) 



where v and q are real numbers (the case of v = is trivial) and q sets the asymptotic position of 
the integration contour. The application of Eq. fl3.1| ) to Feynman integral calculus was initiated 
in Refs. [23, 24] (see also Ref. [28]) and is an algorithmic procedure which can be completely 
automatized, as e.g. in the Ambre.m package [30] in MATHEMATICA. 



In general, the contour in Eq. ( |3.1[ ) is not necessarily a straight line and its standard definition 
is such that the poles of T(v + z) (at z = — i — v with i being non-negative integer) are all to the 
left and the poles of T(—z) (at non-negative integers) are all to the right of it. The condition on 
the poles of the T-functions can be satisfied by such a contour in Eq. fl3.l|) if and only if q < 
and v > 0. However, as a key observation, Ref. [27] realized straight-line contours parallel to the 
imaginary axis in an algorithmic way. If v < 0, we start with a curved contour that fulfills the 
condition on the pole and then deform it into a straight line taking into account the residua of 
the crossed poles according to Cauchy's theorem. This procedure lends itself to implementation 
in computer codes for the evaluation and manipulation of MB integrals, such as in the MB.m 
package [31]. 



-v-2 -v-1 



Im z 




2 Rez 



-v-2 -v-1 



Im z 



2 Re z 



Figure 2: The deformation of a curved contour into the sum of a straight line and a circle for v = —1/2 
and q — —1/4. 
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For instance in Eq. Q3.1| ), if v = —1/2, a possible good choice for q is q = —1/4. The curved 
contour depicted in Fig. Q on the left fulfills the conditions on the poles of the T-functions. Contour 
deformation results in a straight line and a circle around the pole at z = —v = 1/2 as shown 
graphically in Fig. § on the right. Therefore, the MB representation of (a + b) 1 / 2 in terms of a 
vertical straight line contour is given by 

1 y — 1/4+ioo J 

( a + 6 ) 1/2 = fTT7^ / —a 1 / 2 - z b*r(-l/2 + z)T(-z) + Vb, (3.2) 

r(-i/2) y_i /4 -ioo 2vri 

where the first term corresponds to the integral along the straight line with q = —1/4 and the 
second one to the integral along a circle surrounding the pole in z = 1/2 and evaluated according 
to Cauchy's theorem. Note that Eq. ([O]) is not valid for negative integer values of v because of 
r(i/) in the denominator. In these cases we use the binomial expansion 

— — = ( ) a l b~ u ~ l , for v being negative integer. (3.3) 

(« + 6)" few 

We use Eqs. ( |3.1| ) and ( |3.3| ) to convert all sums in the integrands of the soft, collinear and 
iterated integrals of Sect. [2] into products. Then, we apply the relations 

(l-x) = f dyyS(l-x-y), (3.4) 
Jo 

and 

pl n / n \ n / n \ 

/ n^r" 1 *( 1 -E^)=n r (p*)/ r (E^)' ( 3 - 5 ) 

Jo t=i v 3=1 J r=i v i=i / 

to obtain a representation of the original integrals in terms of MB integrals where all the integra- 
tions over a, v, cos(#), cos(0) are performed and only complex integrations along straight lines 
parallel to the imaginary axis are left following the procedure discussed below Eq. fl3.1|) . Upon 
deformation of the curved complex contours all singularities in e are extracted so that it is safe to 
expand in e around zero before doing the complex integration. In this way, the MB representa- 
tions of the required coefficients of the Laurent expansions of the integrals of Sect. [2] are obtained. 
In the next step we convert the complex contour integrations into harmonic sums using Cauchy's 
theorem and finally we evaluate the sums. For the computation of all the harmonic sums we have 
used algorithms for harmonic and nested sums of Refs. [32, 33] as implemented in the XSummer 
package [34]. Typically, symbolic summation of single-scale nested sums leads at intermediate 
stages of the calculation to harmonic polylogarithms (HPLs) (see Ref . [35] for a definition) . How- 
ever, in all cases where analytic result have been obtained by summing series of residues, the 
HPLs could be converted to standard polylogarithms (see also Ref. [21] for a discussion of the 
class of functions appearing in the integrated real- virtual counterterms). 

As an example let us consider the following integral: 

£(x; e, d ) = x 2 [ da a" 1 "' (1 - a) 2d ° [a + (1 - a)x]~ 1_£ [2a + (1 - a)x}' 1 , (3.6) 
Jo 
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which is a typical contribution to the collinear integrals defined in Eq. ( |2.5[ ). The integral in 
Eq. ( |3.6| ) is clearly divergent in the limit e = due to the factor a -1 " 6 in its integrand. The first 
step is to write the MB representation of the integral in Eq. ( |3.6| ). To do this we use Eq. ( |3,1| ) 
twice then Eqs. (|3.4j) and (|3.5|) . We obtain 

ci a ^ r i+io ° dz i r 2+io ° d ^ *2 

/ 2m / ■ 27T? 

xr / —21, -Z2, 2d — 1 — e — 21 — z 2 , 1 + e + zi, 1 + 22) -e + zi + z 2 
X 1 2d - 1 - 2e, 1 + e 



where we have introduced the notation 



a\, a 2 , • • • , a n 
61 , 62 , . . . , fr m 



IX I " ' lll^,)- (3-8) 



For do > 2 we choose (?i = — 1/4 and ^2 = — 1/8 and curved contours such that the real parts of the 
arguments of all T-functions remain positive on them. Note that this implements the requirement 
that the contour separates the left poles from the right ones. In order to use straight-line contours, 
we must add contributions from two residua: the first is due to the residue coming from the pole 
in Z2 = e — zi and then in the resulting one-dimensional MB integral the second is due to the 
residue in z\ = e. Adding these contributions to the starting representation of Eq. (^^), we find 
the MB representation of Eq. (|3.6|) with e close to zero to be given by: 



£{x;e,d ) = x~ 2e T 



-e, 1 + 2e 
1 + e 



+ 91+100 dzi 2 ,„ Z1 x _ 2er I - Zl , -e + Zl ,l + e + z u l + e-z x 
Iqx-ioo 2m \ 1 + e 



Jqi-ioo 2m Jq 2 —i 00 2m 

xr / -zi, -z 2 , 2d - 1 - e - z x - z 2 , 1 + e + z x , 1 + z 2 , -e + z x + z 2 \ , 
\ 2d - 1 - 2e, 1 + e J ' 1 ' ) 

At this point we see that the singularity in e = is isolated in the first term of this equation. 
In particular the pole comes from the factor T(— e) of this term. This shows that the extraction 
of poles comes out in a very convenient way: in practice we have only deformed contours and 
computed residua. As a matter of fact, this is one of the strong points in the application of MB 
methods to phase space integrals, and the straightforward way of extracting infrared poles has 
been already discussed in Refs. [36,37]. 
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As a next step we can perform the expansion around e = and we obtain 

1 /-91+ioo d 

£(x;e,d ) = — + 21og(x)+ / -A 2~ Zl r(l - z u -z x , z x , 1 + zi) 

e Jot —ioo ^ 7r? 



+ r^ +io ° dzi /^+ io ° d£2 2Z2 

Jqi-ioo ^Tri J -ioo ^iri 



X 



-Zl-Z 2 



<72— ioo 



xr 



-zi, -z 2 , 2d — 1 — zi — z 2 , 1 + zi, 1 + z 2 , zi + z 2 
2d - 1 



(3.10) 



The first integral can be easily computed. We close the contour to the right and compute the 
residua coming from the poles enclosed in it at zi = n; n > 0. The residua are given by 
(l/2) n log(2). Thus, multiplying by an overall minus sign due to the clockwise orientation of 
the contour, we find 



I 



91+100 dzi ^2, 

^ 2- 2i r(l - zi, -zi, Zl ,l + zi) = - £(1/2)" log(2) = -2 log(2) . (3.11) 



51— ioo 



n=0 



Next we evaluate the second integral closing both contours to the left. We begin with the integra- 
tion over the variable zi and we have two different T-functions that contribute with poles. The 
first one is T(l + zi) which exhibits poles in zi = —n; n > 1 and the second one is T(zi + z 2 ) 
which contributes with poles in zi = —n — z 2 ; n > 1. Computing these residua, we obtain for the 
second integral 



/•" 1+io ° dzi /" ?2+ioo d£ 2 

/qi-ioo 27Ti Jq 2 —i 00 2m 
X '■'< 1 - ]X dZ2 

2ri 



-Zi, —z 2 , 2d — 1 — Zi — Z 2 , 1 + Zi, 1 + Z 2 , Zi + Z 2 
2d - 1 



1 J ao — 



1 Jq 2 -ioo 



"I) 



n+1 



(n-1)! 
(-*) 



22 „,fl— Z2 



2 22 x 



n, -z 2 , 2d — 1 + n — z 2 , 1 + z 2 , —n + z 2 
2d - 1 



2 Z2 r 



n! 



2d — 1 + n, 1 — n — z 2 , -z 2 , 1 + z 2 , n + z 2 
2d - 1 



. (3.12) 



Now we can do the remaining integration over z 2 . In this case the poles of both the integrands 
are in z 2 = — m; m > 1 and the corresponding residua are given by 



£(*;e,d ) = -- + 21og(|)-log(2) £ (-) W 

m,n=l 

oo / 

- E (!)"> 

m,n=l N 



2d - 2 + n 
n 



(3.13) 



__ n f 2do — 2 + m + n 
m + n 



Si(2do - 2 + m + n) - 5i(m + n) + log (|) 
where the harmonic sums Si(n) are defined as [32,33] 

n 

Si(n) = Y^- = i,(n+l)+ lE , (3.14) 



1=1 
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with 7£ being Euler's constant and tp(x) being the polygamma function, i.e. the first derivative 
of the logarithm of the T-function. The first double sum amounts to 



m,n=l 



1 



1 



(1-x) 



2d -l 



(3.15) 



We are not able to perform the second summation for arbitrary do- However, choosing integer 
values, do > 2, these sums simplify significantly. Indeed, if do is a positive integer, then both 



2c?o — 2 + m + n 
m + n 



2c?o — 2 + m + n 
m + n 



[Si(2d -2 + m + n)- S^m + n)] , (3.16) 



are polynomials in m and n. This implies that the double sums in Eq. ( |3.13| ) can be written in 
terms of the functions 



Li*: (a;) 

(l_a;)l-fc 2^i=0 

where Li/%(x) are the classical polylogarithms [38] and 



n=l 



' k )x- k ~ i 
i I 



if k > , 
if fc < , 

are the Eulerian numbers: 



(3.17) 



i+l 



" ^ r )( i • /,■<(). 



(3.18) 



Therefore, Eq. ( 3.13| ) becomes for example with the choice do = 2, 



1 , ,J 1 A x 2 (3x 2 - 15x + 14) 

£ ' do) = ~~e + bg(2) V 1 " O^W J " 



(x 6 - 9x 5 + 33x 4 - 78x* + 108x 2 - 72x + 16) 



+ 



^l-x) 3 (2-^ 3 



X 



log (-)+ 0(e) 



(3.19) 



Looking at this expression we notice that even if the integral in Eq. ( |3.6[ ) is well defined for x = 1 
or x = 2 some of its individual terms diverge in these limits. Nevertheless, the full result has a 
well defined limit in x = 1 or x = 2. Indeed, 



1 53 

lim£(x-e,d ) = — + — - 161og(2) + 0(e) , 

e 6 



lim £(x; e, do) 

x^2 



1 _ 8 

e ~ 3 



+ 21og(2) + 0(e). 



(3.20) 
(3.21) 



This completes the discussion of our example and demonstrates that the Laurent coefficients 



are given by simple functions in x only. Looking back at how Eq. (3.1) has enabled us to arrive 



at Eq. (3.7) starting from Eq. (p]q) it is obvious that more complicated integrals such as nested 
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ones denned in Sect. 2.2 result in increased numbers of Mellin integrations and shifted arguments 
of the T-functions. Also, in the case of nested integrals the order of the singularities is higher. 
However, the extraction of the poles in e always reduces the dimensionality of the MB integrals 
(in our example from two to zero). Hence, in general if we start from a high-dimensional MB 
representation, the contributions to the poles' coefficients have a much lower dimensionality of 
the Mellin integrals, which allows for an analytic computation of the coefficients of the poles in 
the e expansion even for the most complicated integrals. This example also shows that for the 
analytic computation by means of a MB representation one should choose do,d' to be positive 
integers and transform the regions of integrations in the integrals defined in Sect.|| to [0, 1]. In 
this work we simply choose «o = Ho = 1 an d consider the cases do = d' = 2, 3 which as discussed 
in Ref. [22] are the natural choices for the infrared subtraction for processes with two and three 
outgoing jets respectively. Nevertheless, we stress that in principle any choice of do,d' Q > 2 can 
be used in a computation of m-jet production, for any m. Thus there is no need to recompute 
any integrals even for processes with more than three jets. Furthermore, the appearance of a 
hypergeometric function in the integrand of Eq. ([0]) as happens for example in the last row of 
Tab. |l| does not essentially change the complexity of the computation. The reason is that the 
hypergeometric function 2-F1 has a simple MB representation: 



n i » r +i °° dz ( V* r c,a + z,b + z,-z\ 
,fi(o,6,q«,)= / ^(-«0r b +g > (3-22) 



q— ioc 



where the integration contour separates the poles of the T(- ■ ■ + z) functions from the poles of the 
r(- • • — z) function as usual. 

In closing this Section, we would like to mention another virtue of the MB method. For a 
given phase space integral of Sect.0, the corresponding MB representations show good convergence 
properties if evaluated numerically along the complex contours. Thus, the multidimensional 
numerical integration of MB integrals, such as in Eq. ( |3.9|) is straightforward with the help of the 
CUBA library [39], which provides an independent check. Moreover, it also presents a quick and 
reliable way of obtaining numerical results for the (smooth) 0(e°) terms in the Laurent expansions 
of all integrals for the real-virtual counterterms in the paper. 



4. Collinear integrals X 



In this Section, we show the analytic results for the collinear integrals defined in Eq. (|2.5| ) for 
which the case k = is needed only for the first row in Tab. p] and k = 1 is needed for all of them. 
Analytic expressions for the first two cases of Tab. [I] have already been computed in Ref. [21]. 
Here we fix do = 3 and give the explicit expressions for this case as an illustration of the form of 
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our results. For the Laurent expansion we obtain 



J(x;e;a = 1, d ; k, k, 5, g^) 



S k , 



2(2 - 6) e 2 

:(±) 



2S k ilog(x) 



+ 



l-S k , 



3-5 2[l + fc(l-<y fc ,_ 1 )] 



+ S KA g^>{x) + F(x; e, d ,k) + 0(e), 



(4.1) 



where I = A,B,C,D (see Tab. [T]), k = — 1,0, 1,2 and 8ij is the usual Kronecker 5. Here we have 
introduced the two functions Qj k (x) and Tix; e, do,k). The function Q^fj} is a matrix in / (rows) 
and k (columns) defined as follows: 



(±), 



x 



( 


§C 2 + ±log 2 (x) 


1 


1 

2 




13 \ 
36 




±f)C 2 +(±T±)log 2 (x) 


1 




1 
1 


13 _|_ 1 
16 4 




(|±i)C2 + ilog 2 (x) 


1±I 




3 
8 


13 _i_ 11 
36 31 36 


i CM 

\ 1 36 


T^)C2+(^±^)log 2 (x) 


l±i 




1 

8 


36 ^ 18 / 



(4.2) 



and choosing e.g. do = 3 for the function F{x; e, do, — 1), we obtain 



.F(x;e,d = 3,-1) 





e,do 


= 3,0) 




e,do 


= 3,1) 


J^(x- 


e,d 


= 3,2) 



O 1 

--C 2 + log 2 (x) - ^^Sf-lC^; 35 > " 133 > 188 > " 116 > °. °) 



~P {5 \x; 49, -193, 281, -173, 24,0) + P^ir: 1. -5. 10. -10.5.2) . ( I. I.) 



,(5) 
2,- 



(4.3) 



-r^^f-i^; 25 > - 116 > 212 > ~ 192 > 96 > °) - p 2,-i( x ; !» - 5 > 10 > 5 > 2 ) 



i 



jF(x; e, do = 3, 0) 



(4.5) 



80x 



log(2) + 



(1 -x) e 



Pffifa 51, -861,6523, -29212, 



3(2 -x) 6 ™ GV ~ y ' 36(2 -x) 6 ' 

85505, -171607, 241761, -240096, 164864, -74000, 19120, -2048) 

+ iro~ X \e p i^ ( X] x > ~ 17 ' 130 ' ~ 590 ' 1765 ' ~ 3734 ' 5748 > -6360, 



4880,-2480,784,-128). 
Here we introduced the short-hand notation 



(4.6) 



M) (rr . n {k) „(*)<> _ Li n (l-x) sr 



r n,k V X ' a m i ■ ■ ■ i u > 



(1 - x) r 



i=0 



(4.7) 



According to their definition, the limit of the functions given in Eqs. (|4.3[)-(p~6|) must be finite in 
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x = 1 even if some terms are separately divergent. Indeed computing the limit at x = 1 we find 



lim F(x; e, do = 3, — 1) 



x-*l 



lim F{x; e, do = 3, 0) 

X— >i 

lim e, do = 3, 1) 
lim e, do = 3, 2) 

X— >1 



8731 

3600 
257 

~60" ' 
257 

120 ' 
1801 

90 



<2, 



+ ?log(2) 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



In Fig. || we compare the analytic and numeric results for the e° coefficient in the expansion 
of I(x, e; 1, 3; 1, — 1, 0, g^) and I(x, e; 1, 3; 1, — 1, l,gffl) for k = —1, ao = 1 and do = 3 as 
representative examples. The agreement between the two computations is excellent for the whole 
x-range. The numeric results have been obtained using standard residuum subtraction and a 
Monte Carlo integration program as explained in detail in Refs. [21,22]. This shows that the 
expansion coefficients of all the collinear integrals X and hence also of the collinear subtraction 
terms are smooth functions of the kinematical variable x. 



o 4.0-10- - 



2 10 
I 



C-type, k = 1, k = —1, 8 = 




• ao = 1 



□ 1-cr □ 2-cr □ 3-tr 



-8 -7 -6 -5 -4 -3 -2 -1 



D-type, k = 1, k = — 1,6 = 1 




♦ a = 1 



□ l-CT □ 2-tr □ 3-tT 



-8 -7 -6 -5 -4 -3 -2 -1 



Figure 3: Representative results for the C-type and D-type integrals. The plots show the coefficient 
of the O(e ) term for k = -1 in l(x, e; 1, 3; 1, -1, 0, g^) (left) and l(x, e; 1, 3; 1, -1, 1, g^) (right) with 
do = 3 and ao = 1. 

The complete results for all necessary cases (like in the later Sections) are of considerable 
size, such that we shall not list them here. They are all contained in a MATHEMATICA file provided 
with the sources of the paper on the archive http : / / arXiv . org. 



- 15 - 



5. Nested collinear-type Z*Z and X*J integrals 



In this Section we discuss the analytic computation of the nested collinear integrals defined in 
Eqs. (plOD-CT- 

As an example we show explicitly the fully analytic result for the case X*X r (x, e; 1, 3; —1, 2) 
for which we were able to compute the complete pole structure analytically Choosing do = 3 and 
«o = 1 we get 



X*X r (x, e; 1, 3; —1, 2) 



11 ( 2 1 . , 
"l2^H"9 + 3 l0g(X ' 

+~ log(l - x) log(x) + ~ U 2 {x) \ + (1 _^, /2)5 



1 



+ 



1 



(1 



1 



c 2 



25 
36 



log(x) 



+ 



1/6 



+- 



(l-x/2) 4 
1 



+ 



1 



(1 - xf 



1 



log(x 

72 18 &v ' 

1/18 



+ 



— log 

6 6 

1/12 



x 

r ) I + 



(1 



13 1 , / \ 
"36 + 6 l0g(x) 



1 2 , i 

log(x) + / ■ ■> 

6 9 6V V (l-x/2) 2 



(l-x/2) 3 
1 



(1 



25 7 i ^ 
log X 

72 12 sv ' 



1/24 



31 1 19 2 2 

+- — ! — — H 1 — log(2) H log(x) H — log(l — x) log(x) log (x) 

(l-x/2) 216 6 w 9 v ; 3 V ; V ; 3 5 v ' 



+ 3^2^) 



- + 0(e°). 



(5.1) 



This result is representative, because its form is typical of all the collinear nested integrals. The 
plot of the 0(e _1 ) coefficient of this Laurent expansion for the integral 2*T r (x, e; 1, 3; —1, 2) is shown 
on the right side of Fig. ||] together with the comparison with the numerical evaluation obtained 
using standard residuum subtraction and Monte Carlo numerical integration. On the left side of 
Fig. U we plot the same coefficient of the Laurent expansion for the integral X*Xj(x, e; 1, 3; —1, 2). 
For both cases we note that the agreement between the numerical evaluation and the analytic 
result is excellent. These plots show also that the coefficients of the Laurent expansion of the 
nested collinear integrals 1*1 are very smooth functions of x. 



We note that in the Laurent expansion of I*I r (x, e; 1, 3; — 1, 2) in Eq. (5.1) there are some 
terms that are divergent in x = 1. However according to its definition in Eq. ( |2.11| ) the limit in 
x = 1 must be finite. To verify this is a further check of the correctness of the result. For the 
case of Eq. (5.1) we obtain that: 



limM r ( I , e ;l,3;-l ) 2)=-ii-: j + 



2 1 

7 2 



3091 
675 



31 



+ ^2-y log(2)) - + 0(6 U ). (5.2) 



The case of X*X r (x, e; 1, 3; 2, —1) is more difficult. For this integral we are unable to compute 
the coefficients of the e poles in a fully analytic form. The reason is that in its Mellin-Barnes 
representation also three-fold MB integrals are involved. For this case the coefficient 0(e~ 3 ) and 
0(e~ 2 ) are fully analytic but the coefficient of 0(e _1 ) is semi-analytic. This last coefficient is 
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k = — 1, I = 2 



-5.0- lO- 
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Figure 4: Representative results for the Z*Z-type integrals. The plots show the coefficient of the 0(e 1 ) 
term for k = — 1 and I = 2 in T*Ti(x, e; 1, 3; — 1, 2) (left) and I*I r (x, e; 1, 3; — 1, 2) (right) with do = 3 and 
a = 1. 



thus written in terms of an analytic expression to which a three-fold MB integral must be added. 
The remaining MB integral can be efficiently computed in MATHEMATICA by use of the package 
MB.m [31]. Explicitely for I*I r (x, e; 1, 3; 2, —1) we have: 



l*I r (x,e; 1,3; 2,-1) = -\ \ + 



(1 - x/2) 6 



12 l0g (I 



+ 



1 



(1 - x/2) 5 
5/72 



5 5 , fx 

I log - 

12 12 B V2 



+ 



1/6 



+ 



5/24 



(1-x) 
1/12 



5 l a l °g( X ) 



1/18 5/144 



(1-x) 4 (l-x/2) 4 (1-x) 3 
1/48 



1/24 
+ tt 1 r + 



59 1 M 

(l-x/2) 3 ' (1-x) 2 ' (l-x/2) 2 ' (1-x) ' (l-x/2) "72 + 2 ° g[X) 
1 /25 21 5 5 21 

(1 _ x/2)6 log(2) + 4 log2(2) + 3 log(2) log(1 " X,2) + ~% l ° g{x) 



5 5 5 5 

~2 lo s( 2 ) lo s( x ) + — log(l - x) log(x) - - log(l - x/2) log(x) + - log 2 (x) 



l Li2 (f) + il Li ^) + (r^F 



\ C 2 - \ log 2 (2) - \ log(2) log(l - x/2) 
3o 3 



— log(x) + i log(2) log(x) + i log(l - x) log(x) + i log(l - x/2) log(x) 



I log 2 (x) + i L i 2 (I) 4- I Li2 (x)) + Ir --L^ 5 ^ 



23 25 , 71 , , . 

( 2 + - log(2) 



-- log 2 (2) - - lo g (2) logtl - ,-/2) - lo g (,-) + - loo(2) log(..-) + 



5 
3 



17 
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-j2 ~ x ) lo s( x ) + | lo § ( X ~ |) lo §( x ) ~ | lo g 2 (z) + | Li 2 (|) - ^ Li2(jc) 

+ (T3^)4 (I ~ \ ^ + \ log ( 2 ) ~ Y log ^ + 1 log ^ ~ x ^ log ^ + 1 Lia ^ 

1 / 59 1 /x\\ 1 (I 1 7 

+ (l-x/2)4 (,-48 + 6 l0g UJJ + (l^^-i6-12 C2 -36 l0g(:c) 

1 1 \ 1 / 29 1 4 

H logfl — x) log(x) H Li 2 (x) + — t-t I log(2) H — log(x) 

12 &v ' ^ ' 12 ^ 7 (l-x/2) 3 V 432 6 w 9 v ; 

1 (I 2 , . , 23 , . A 1 [211 2, /n . 7, 

+ 7 rrr h - log(2) log(x) + -. — rrr log(2) + - fog(x) 

(1-x) 2 V 27 9 5W 36 6V 7 (l-x/2) 2 V864 9 BK J 9 5V v 

1 / 31 59 , , A 1 /139 1. ... 4. . A 1177 7 A 

+7 r log(x) + — log(2) + - log(x) h - C 2 

(1-x) V 72 24 BV 7 (l-x/2) V288 3 w 3 v 7 4 32 8 

- i log(2) + ^ log 2 (2) + i log(2) log (l - |) + ^ log(x) - i log(2) log(x) 

I 1 / x \ 5 I / x \ 1 

+- log(l - x) log(x) - - log (l - -J log(x) - - log 2 (x)) - - Li 2 ^-J + - Li 2 (x) 



+MBint[xl 



^ + 0(e°), (5.3) 



where MBintfx] is a three- fold Mellin-Barnes integral, which for this case is given by 



MBintfxl 



gi — ioo 27TZ Jq 2 -i OQ 2ixi Jq 3 _ioo 27ri 



p / -Zl, 1 + «l, 3 - «2, -2 + z 2 , 5 - zi - z 2 - z 3 , -z 3 , 2 + z 3 , z x + z 2 + z 3 \ 
{ 4, A-z 2 J ' 



where q± = q 2 = q 3 = — 1/4. 



Similarly to the analytic expression of Eq. (5.1), also in this case we have many terms that 
are singular in x = 1 even though the full expression is well defined. Moreover in cases like this 
where we have a semi-analytic expression we find that the analytic part and the remaining part 
expressed in terms of a three- fold MB integral are separately well defined in x = 1. In particular 
for the case of the integral Z*Z r (x, e; 1, 3; 2, —1) in Eq. ( |5.3; ) we obtain the following limit: 

,. rr t / , „ „ \ 1 1 / 607 40 , . A 1 Z77349 509 A 

limZ*X r (x,e; 1,3; 2,-1) = + h — log(2) ) ^ + h C 2 

x-+i rv ' ' ' ' ' ; 6 e 3 V 60 3 7 e 2 Y 144 00 24 s 

:Ul 1 log(2) + log 2 (2) + MBint[lA - + 0(e°), (5.5) 



45 ov ' 3 ° v ' L 7 e 

where MBintfl] is given by 

MBintfl] = 0.329808. (5.6) 



This number is the result of the MB integral in Eq. ( |5.4[ ) with the choice x = 1 obtained using the 
MATHEMATICA package MB.m [31]. Finally we note that this example is representative for a small 
subset of the collinear nested integrals which have these features. They are X*2j(x, e; 1, 3; k, I) and 
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1*tL r (x, e; 1, 3; k, I) with k = — 1, 1, 2 and I = — 1 and T*J (x, e; 1, 3, 1, 3; k) with k = — 1. The results 
for the pole structure of all the remaining cases of nested collinear integrals are fully analytic. 



O 5 - ' 



X*J, k = 




• aa = yo = 1 



□ 1-cr □ 2-cr □ 3-cr 
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Figure 5: Representative results for the I<7"-type integrals. The plots show the coefficient of the 0(e 1 ) 
term for k — in I*3f(x, e; 1, 3, 1, 3; 0) with do = d Q = 3 and ao = yo = 1. 

In Tab. || we list numerical values for the non-trivial coefficients of the e-poles (i.e. the 
0(e~ 2 ) and 0(e~ 1 ) coefficients) of the nested collinear integrals T*T r (x, e; 1, 3; — 1, 2) and T* 
I r (x, e; 1, 3; 2, — 1). These numbers have been obtained using the fully analytic expression in 
Eq. flS.lj ) and the semi-analytic one in Eq. (|3l|). Numbers for the 0(e°) coefficient for the same 
representative integrals are listed in Tab. [3|. In this case they have been entirely obtained evalu- 
ating their MB representations. 

Finally in Fig. || we plot as a further example the fully analytic result for the first order e-pole 
for T* l 7(x, e; 1, 3; 0) together with the numbers obtained numerically using standard residuum 
subtraction and Monte Carlo numerical integration. As for all other cases the agreement is 
excellent and the coefficient is given by a very smooth function of x. 



6. Nested soft-type J*J integrals 



In this Section we discuss the analytic computation of the integrals defined in Eqs. ( |2.13| )-( [2.15| ). 
For them we were able to compute a fully analytic result for the coefficient of the Laurent expan- 
sion up to 0(e~ 2 ). The 0(e _1 ) coefficient is computed semi-analytically similarly to the nested 
collinear integral T*l r (x, e; 1, 3; 2, — 1) discussed in Sect. || As a representative example we show 
the structure of the fully analytic part of the result for the nested soft integrals J*J. For example 
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choosing d' Q = 3 we have: 

1/22 \ 1 1 

J*J lk (Y; e; 1, 3) = -j + I — - 2 \og(Y) \ -3 + W(y) ^ + O^ 1 ) , (6.1) 

J*J ir {Y- e; 1, 3) = i i + _ iog(y)) !+(§-§ log(Y) + log 2 (y) 

+^Li 2 (l-y))l + 0(6- 1 ) (6.2) 

and finally 

1/22 \ 1 / 1 \1 

J*J kr (Y; e; 1, 3) = ? + ( — - 2 log(y) J - + ( W(y) - C2 + - Li 2 (l - F) J ^ + °( e_1 ) • ( 6 - 3 ) 

The function TOY) which appears in Eqs. (|6.1|) and (|6.3[) is given by 



H(V) = £ - 2 6 + ^ + 33^-U^^-44 [og(y) + 2 log2(y) +4L . (1 _ n 

(6.4) 

Also for this function even if some terms are singular at Y = 1, we still have that the limit is well 
defined. Indeed we find 

97 

hmW(y) = y -2C 2 . (6.5) 

For these three soft-type integrals the 0(e -2 ) coefficient has been plotted in Fig.|6| using its 
fully analytic expression Eqs. (|6.l| ) and fl6.4|) and its numerical evaluation obtained using residuum 
subtraction and Monte Carlo integration. The agreement is excellent and the analytic result 
confirms that also the coefficients of the Laurent expansion for the J* J integrals are smooth 
functions of Y . 

The numbers in Tab. ||| have been obtained evaluating the nested soft integral ^*J7ik{Y; e; 1, 3) 
using the fully analytic expression in Eq. fl6.1|) for the 0(e~ 3 ) and 0(e -2 ) coefficients. For the 
0(e -1 ) coefficient a semi-analytic expression in terms of a MB integral has been used and finally 
the representation only in terms of MB integrals has been evaluated for the 0(e°) coefficient. 



7. Nested soft-collinear K,*J integral 



In this last Section we discuss the pole structure of the integral defined in Eq. ( |2.16 ). In this case 



the result for the Laurent expansion is very simple because the integral has no dependence on the 
kinematics. The coefficients of the poles in lC*J{e, 1, 3) with d' Q = 3 and yo = 1 read: 

^ , , 1 1 11 1 557 1 / 10825 5 \ 1 n . n , . . 

This completes our discussion of the analytic computation of the fundamental integrals that 
contribute to the singly- unresolved counterterms. 
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Figure 6: Representative results for the J^Jf-type integrals. The plots show the coefficient of the 0(e 2 ) 
term in J*Jik(Y, e; 1, 3) (left), J*tfi r {Y, e; 1, 3) (right) and J%Jkr(Y, e; 1, 3) (bottom) with d' = 3 and yo = 1. 



8. Conclusions 



In this work we have completed the evaluation of all integrals needed for the computation of the in- 
tegrated real-virtual counterterms of the subtraction scheme for NNLO jet cross sections proposed 
in Refs. [17-19]. We have discussed representative examples for all types of soft and collinear as 
well as nested integrals in Sects. |I|-|7] (the complete results are contained in a MATHEMATICA file). 
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These integrals (i.e. their Laurent expansions in e to sufficient depth) have to be computed once 
and for all and their knowledge is necessary in order to make the subtraction scheme an effective 
tool. We have achieved this task by deriving MB representations for all integrals under consid- 
eration and, in a subsequent step, we have performed analytically the summation of the nested 
sums over the series of residues. In some cases, this second step of summing the series has not 
been achieved and we have resorted to a numerical evaluation of the MB integrals in the complex 
plane. As a further check, all MB representations for both the numerical and, if available, the 
analytic results have also been compared against an independent evaluation of the integrals using 
standard residuum subtraction together with the Monte Carlo integration in Ref. [22]. We have 
shown, that all integrals contributing to the real-virtual counterterms are smooth functions. For 
practical applications, this means that all integrals (in particular the finite in e contributions) can 
be used in terms of interpolating tables, which are computed once and for all. 

Files of our results can be obtained from the preprint server http : / / arXiv . org by down- 
loading the source. They are also available at [40] or from the authors upon request. 
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LO B10\ X ) 


I*I r (x, e; 1, 3; —1, 2) 


l*l r (x,e; 1,3; 2,-1) 




0(e- 2 ) an. 


0(e- 1 ) an. 


0(e- 2 ) an. 


0(e 1 ) semi-an. 


-10. 


-7.89751 


-374.957 


-15.9061 


-759.736 


-9.66667 


-7.64166 


-351.104 


-15.3944 


-711.688 


-9.33333 


-7.38582 


-328.036 


-14.8828 


-665.211 


-9. 


-7.12998 


-305.753 


-14.3711 


-620.305 


-8.66667 


-6.87413 


-284.256 


-13.8594 


-576.969 


-8.33333 


-6.61829 


-263.544 


-13.3477 


-535.205 


-8. 


-6.36245 


-243.618 


-12.836 


-495.012 


-7.66667 


-6.10661 


-224.477 


-12.3243 


-456.389 


-7.33333 


-5.85076 


-206.122 


-11.8126 


-419.337 


-7. 


-5.59492 


-188.552 


-11.301 


-383.857 


-6.66667 


-5.33908 


-171.768 


-10.7893 


-349.947 


-6.33333 


-5.08324 


-155.769 


-10.2776 


-317.608 


-6. 


-4.82739 


-140.556 


-9.7659 


-286.841 


-5.66667 


-4.57155 


-126.128 


-9.25423 


-257.644 


-5.33333 


-4.31571 


-112.485 


-8.74256 


-230.019 


-5. 


-4.05986 


-99.628 


-8.2309 


-203.965 


-4.66667 


-3.80402 


-87.556 


-7.71928 


-179.482 


-4.33333 


-3.54818 


-76.269 


-7.20772 


-156.573 


-4. 


-3.29234 


-65.7665 


-6.69628 


-135.237 


-3.66667 


-3.03649 


-56.0477 


-6.18508 


-115.476 


-3.33333 


-2.78065 


-47.1111 


-5.6743 


-97.2944 


-3. 


-2.52481 


-38.9536 


-5.16432 


-80.6957 


-2.66667 


-2.26896 


-31.5702 


-4.65576 


-65.6862 


-2.66666 


-2.01312 


-24.9522 


-4.14969 


-52.2723 


-2. 


-1.75728 


-19.0853 


-3.64776 


-40.4585 


-1.66667 


-1.50144 


-13.9478 


-3.15236 


-30.2408 


-1.33333 


-1.24559 


-9.50936 


-2.66658 


-21.5965 


-1. 


-0.989751 


-5.73082 


-2.19382 


-14.4712 


-0.66667 


-0.733908 


-2.5675 


-1.73699 


-8.76877 


-0.33333 


-0.478065 


0.0265877 


-1.2975 


-4.35339 


0. 


-0.222222 


2.09462 


-0.874704 


-1.06702 



Table 2: Numerical values for the 0(e -2 ) and 0(e _1 ) coefficients of 1&E r (x, e; 1, 3; — 1, 2) (second and third 
column) and T*T r (x, e; 1, 3; 2, — 1) (last two columns) for various values of log 10 (x) (first column). These 
numbers have been obtained evaluating the fully analytic expression in Eq. (5.1) and the semi-analytic one 
in Eq. dJ) 
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togiofc) 


T*T r (x,e;l,3;-l,2) 


l*l r (x,e; 1,3; 2,-1) 


-5. 


-1643.45 


-3380.25 


-4.66667 


-1354.81 


-2792.08 


-4.33333 


-1104.32 


-2276.69 


-4. 


-886.741 


-1829.25 


-3.66667 


-699.713 


-1444.98 


-3.33333 


-541.331 


-1119.05 


-3. 


-409.041 


-846.661 


-2.66667 


-300.305 


-622.985 


-2.33333 


-212.59 


-443.178 


-2. 


-143.341 


-302.315 


-1.66667 


-89.9699 


-195.384 


-1.33333 


-49.9194 


-117.263 


-1. 


-20.7583 


-62.7773 


-0.66667 


-0.267788 


-26.8566 


-0.33333 


13.4889 


-4.81253 


0. 


22.1523 


7.37736 



Table 3: Numerical values for the 0(e°) coefficient of T*T r (x, e; 1, 3; — 1, 2) (second column) and X* 
T r (x, e; 1, 3; 2, — 1) (last column) for various values of log 10 (x) (first column). These numbers have been 
obtained evaluating their MB representation. 
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logioon 


J*J ik (Y;e; 1,3) 




0(e" 3 ) an. 


0(e- 2 ) an. 


0(e 1 ) semi-an. 


0(e°) MB 


-10. 


53.385 


1430.99 


25680. 


347094. 


-9.66667 


51.85 


1350.22 


23545.6 


309328. 


-9.33333 


50.3149 


1271.81 


21533.4 


274744. 


-9. 


48.7799 


1195.75 


19639.8 


243157. 


-8.66667 


47.2448 


1122.05 


17861.1 


214389. 


-8.33333 


45.7098 


1050.7 


16193.8 


188265. 


-8. 


44.1747 


981.714 


14634.2 


164617. 


-7.66667 


42.6396 


915.081 


13178.6 


143283. 


-7.33333 


41.1046 


850.805 


11823.6 


124106. 


-7. 


39.5695 


788.886 


10565.3 


106934. 


-6.66667 


38.0345 


729.322 


9400.38 


91621. 


-6.33333 


36.4994 


672.116 


8325.04 


78027.5 


-6. 


34.9644 


617.265 


7335.7 


66018.2 


-5.66667 


33.4293 


564.771 


6428.75 


55463.9 


-5.33333 


31.8942 


514.633 


5600.58 


46240.8 


-5. 


30.3592 


466.852 


4847.56 


38230.9 


-4.66667 


28.8241 


421.427 


4166.08 


31321.6 


-4.33333 


27.2891 


378.358 


3552.51 


25405.6 


-4. 


25.754 


337.645 


3003.24 


20381.7 


-3.66667 


24.219 


299.287 


2514.63 


16153.6 


-3.33333 


22.6839 


263.283 


2083.06 


12631. 


-3. 


21.1488 


229.632 


1704.89 


9728.88 


-2.66667 


19.6138 


198.33 


1376.45 


7367.63 


-2.33333 


18.0787 


169.37 


1094.05 


5473.13 


-2. 


16.5437 


142.739 


853.961 


3976.65 


-1.66667 


15.0086 


118.417 


652.369 


2814.76 


-1.33333 


13.4736 


96.3641 


485.392 


1929.49 


-1. 


11.9385 


76.5204 


349.046 


1268.34 


-0.66667 


10.4034 


58.7892 


239.262 


784.581 


-0.33333 


8.86839 


43.0286 


151.932 


437.509 


0. 


7.33333 


29.0435 


82.998 


192.684 



Table 4: Numerical values for the 0(e" 3 ), 0(e" 2 ), 0(e" 1 ) and 0(e°) coefficients of J*J ik (x; e; 1, 3) for 
various values of log 10 (Y"). The numbers have been obtained from Eq. ( |6.l[ ), the semi-analytic one for the 
0(e _1 ) coefficient and MB integrals for the O(e ) coefficient. 
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